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ABSTRACT 


Wieeins, Edwin George, M.S. in Nuclear Engineering, Purdue 


University, January 1969. An Investigation of the Use of the Hot-Wire 


Anemometer in Non-Isothermal Air Flow. Major Professor: Alexander 
sesonske. 


Procedures were developed for the calibration of the hot-wire 
anemometer in a non-isothermal flow. Various methods of using the re- 
sults of calibration to compute turbulent fluctuations of velocity and 
temperature were considered, and the Arya and Plate modification of the 
Kovasznay Fluctuation Diagram method was found to give the most accurate 
results. The method requires independent meaSurement of the temperature 
fluctuations which was obtained with a rapid response thermocouple. 

The calculated values of velocity turbulence intensity show devi- 
ations of 36 percent from the results of Laufer; however, this deviation 
Beeeo be expected in view of the estimated experimental error. These 
errors were probably caused by contamination of the wire during the 
Menegthy calibration process. 

All measurements were made in air flowing in a 1.0625 inch diameter 
pipe at Reynolds Numbers between 30,000 and 50,000. A heat input of 430 
watts over the last 58 inches before the probe was used to generate 
temperature differences between wall and centerline of 48 to 63°F at the 
probe location. 

The procedures developed and the experimental problems identified 
are the contribution of this investigation to a comprehensive research 
program of non-isothermal hot-wire measurements in fluids of a wide range 


of Prandtl Numbers. 














INTRODUCT LON 
The hot-wire anemometer has been uSed extenSively for the measure- 
ment of velocity fluctuations in isothermal flow of gases. The principles 


? 


: 3 
and procedures involved are discussed by Kovasznay and by Kronauer 
. 4 oe. 

Corrsin has shown that in addition, the hot-wire anemometer has the 
potential to measure Simultaneous velocity and temperature fluctuations 
in non-isothermal flow. Experimental work of this nature has been done 

; 2, ve / 8 
by Corrsin and Uberoi , Gibson Chen and Lin , Deissler , Kunstman , and 
9 
Arya and Plate . The present work was concerned with the development of 
experimental procedures and the identification of problems associated with 
making of such measurements in a non-isothermal air flow. 

In isothermal flow the only flow variable affecting the anemometer 
output is the velocity. Thus the anemometer voltage fluctuation is 
directly proportional to the velocity fluctuation and the anemometer 
steady voltage is directly proportional to the mean velocity. 

In non-isothermal flow, however, both the velocity and ambient 
temperature fluctuations affect the anemometer output. Considerable dif- 
ficulty has been encountered by previous workers in separating the two 

8 ; 
effects. Errors up to 750 percent’ have been reported in non-isothermal 
measurements, while in isothermal measurements 10 percent error or less 
Vs, expected. 
The non-isothermal flow measurements of the present work were pre- 


ceeded by velocity fluctuation measurements in isothermal flow with a 





hot-wire paralleling the work of eee and Sencboun” and temperature 
fluctuation meaSurements in non-isothermal flow with a rapid response 

thermocouple. Thus the velocity and temperature fluctuation intensities 
measured with the hot-wire in non-isothermal flow could be compared with 


Besuits of more direct measurements. 





THEORETICAL BACKGROUND 


The constant temperature hot-wire anemometer consists of a fine 
tungsten or platinum wire maintained by a feedback system at some constant 
temperature above that of the fluid in which it is inserted. The con- 
vective heat loss from the wire is balanced by joulean heating of the 
wire by the current Supplied by the feedback system. 

12 

King solved the convection problem by assuming potential flow 

around the wire and obtained the following relationship: 


7 T -T. 
Nu = (A + BRe*) (1+ 5 ~~ 


where Re = Dup 
rv 


Although it is recognized that King's assumptions are unacceptable, the 
functional form of his relationship has proved to be a satisfactory re- 
presentation of experimental data. Thus it iS Seen that the joulean 
heating current (or voltage) required to maintain the wire at a constant 
temperature is a function of the instantaneous velocity of the flow and 
mreedifference in temperature between wire and fluid. 

In isothermal flow, the temperature difference is constant and the 
System can be used to measure the mean and fluctuating components of the 
instantaneous velocity. In non-isothermal flow the fluid temperature 
has a fluctuating component that affects the heat transfer in addition 


to the effect of the fluctuating component of the velocity. Thus in 





meenciple, tie system can be used to measure velocity and temperature 


fluctuations, since it is sensitive to both fluctuations. 


Hot-Wire in Isothermal Flow 
ee 
The following parales <Gortsm), > treatment for a constant current 


mOL-=wire anemometer. For heat transfer to flow over a cylinder: 





Ie AN ™ 
Nu = (A + B en a) +58 es (1) 
a 
In the range of interest, 
io | 
_ ee a (2) 
2 a 
a 
Whus Equation 1 may be simplified to 
Nu =A+B poe (3) 
By detinition, 
ep, a | mo) 
ek kT = tT) 
W a 
mor CGlectrical resistance heating, 
gi! = ae 
R mDR (5) 


substitution of Equations 4 and 5 into 3 yields 


i 


Ee } AG 1/2 
Ga) eA) + (akan) Re om 
Ww WwW a 


By definition, 


Ee = oe C7) 
vi 





Thus, 


Be 


R(T -T_) 
wow a 


eee 2 
= (mk2A) + (rkxB [ = ] ye 
z (8) 


The quantities in parentheses on the right side of Equation 8 are 


Ponstants. For 


isothermal flow 





constant temperature operation of the hot-wire in an 


the following mean and deviation quantities are used: 





u=udtu! and E=E+e', 
ee) 
Substitution of Equations 9 into Equation 8 yields, after small 
terms have been discarded; 
=) = ie? 12 
Peet 2Ee' _ ed el FO fee 
R(T -T) > (mk&A) + CokQB [| . | yu (1 + = ) 
ow «A 
(10) 
By definition, 
ie - (R_ - RD). 
sas (11) 
According to the binomial theorem 
12 
oy ce 
- = IS 
ra =) a a (12) 
' 
bince = <= |, it 1S assumed Chat higher order terms may be 
meplected in Equation 12. Thus Equation 10 becomes 
Wi a 2 
iE 2Ee' Tks ee oo melee u! 
ee? ano a ce 8) ONG +). 
R = i ( Es 
BCR Rk) R CR R» aR aR ntl Crs) 
Time averaging of Equation 13 yields 
awd | iy? 
E kg k& : = 
Se ae ya (14) 





irk -R) (oR 
W W a O 





aR U 
O 





Subtraction of Equation 14 from 13 gives 








Se 1 ey ' 
Zine ai Tkk pp 22 ] ) u . 
R(R -R) ane Ls: aytl? Ci) 
Ww oOW al 
or 
__ Ee 
—1/2 eee (16) 
R CR Ne Cu u 
Note that 
Bexke , , od 
aR | 7 : Glee 


aot 
end that this same quantity appears as the coefficient of u* in 


Equation 14. Squaring, time averaging, and taking the square 


meor Of Equation 16 yields 





4EY aie = "A ue 
awe vy —= 
fe -R) co’ u oe 
wow a 
Rheplcencoponse Thermocouple in Non-itsothermal F low 
The behavior of a thermocouple over a fairly large range of 
temperature may be expressed as 
E = ES ea 4 (19) 
The following fluctuating quantities are introduced: 
E=E+e' and T=T+t'. (20) 


Sm@pstitution of Equations ZOPinto 1? yields 


PE ese’ = Ey + a(T + ee ee C2) 





Time averaging of Equation 21 yields 


PSE al (22) 


miibtracting Equation 22 from 21 yields 


cy 
e' = at! 
Squaring, time averaging, and takine the square root of Equation 23 
eS = ——_ nea 
"4 ap Seay t'- ‘ a) 


Hot-Wire in Non-Isothermal Flow 
; ; 4 
The following parallels Gorresan > ocbeatment for a constant 
current hot wire anemometer. 


Equation 8 is still valid, however, OP and k are no longer constants. 


Bebstitution of Equation 11 into 8 yields 





2 1 
ee ey, tA gy 1/2, 1/2, 1/2, 1/2 
R (R-R) ‘aR aR “Cu Coes p (2°) 
W Ww a O O Dp 
It is assumed that 
eee 
Gal iG 


Pp 
is constant, and that os teseli ots a1 so constant. 


The following mean and deviation quantities are considered: 


a a a ; (26) 





Substitution of Equations 26 into 25 yields 


=. 2 = 
Ere (ay (k tk’) + 
Mee -R-r') © 0 
teow aia 1/2 € 1/2 27) 
+ (SR Py oy toy A uty ary? 
O 


According to the Binomial Theorem, 


oe aie 1 Gee a 
Peer - Rk) ° =a = i ee eee a 
WoW 31 zl R (R -R ) R -R 
Woo a 





(28) 





Note that 


Wee) > = Eo ne! . Co 


Using Equations 28 and 29 for the left side of 27, and expanding 


the right side of 27, and neglecting small quantities yields: 


E Be 2Ee' 12 
ReCR -R ) R (R -R)* R (R en: 
WwW WwW a Ww WwW a W WwW a e 





A) (k +k!) + (30) 


ee ee ly 2 
g —— —— ae — ——— 
0 SS B) (Qpuk + puk!' + pu'k + @tu ye 
O 


But 


=a oe a = pane ee Ae 1 ! 1 . 
(ouk + puk’ + ou'k + ata! @ x io oe ale= + a pa) (31) 
}/ 














Thus Equation 30 becomes 
=e 
=? ie o = 
a 2-e' ae 
—~. + a = (> A)(kK + k') + (325) 
R (R -R_) R (R -R_) R (R -R ) oO 
mw a Ww a w WwW a 
i 2 Goh) 2 
T2D _p ——, 1/2 sere us (ome 
+ —— oa —_— a 
CR pl DOW At ZR tty YD. 
The time average form of Equation 32 is 
=? 1/2 C 2 
E - 12 = 12D ie) —— 1/2 
= CR ee aR PE B) (puk) ° (33) 
R (R -R_) O O 
ww a 
Subtraction of Equation 33 from 32 gives 
a a 2Ee' mg 
ce > tC Ak (34) 
R (R -R ) R (R -R ) O 
ow (a W ow a 
] ee 
a? oo, gy Couky !? ( a ul i pt ) 
aR, Pr 2 k u 0 ; 
The following relationships are known: 
; ' eet 
= ! ! _ t 
- = Ro i k = nt e = =- 0 T (a5) 
a 
Thus Equation 34 becomes 
EC AA oF i 
S 
= 5 ee C!S, ara) e - 
R AR HR.) RR -R,) O (36) 
ij 2 12 ——— 
, ct ee ee thy U5 
aR Pr 2 k ib = : 
u 


, Equation 36 becomes 





10 


ppl /2p 1/2pet/2, 1/2 


2E | e' = : Pp u 
F (R_-R_) an Pr a ] 
w WwW oO 





i = - en us 
ee R (R -R o (37) 
Ww wW a 


Equation 3/7 may be written 


Ree eps Cay pn 
ww a p 


e' = u! + 
meee Pea 
O 
am R (RO -R) norwlaA -ER a | 
iz O | 
wa 


The coefficient of u' represents the sensitivity to velocity fluctuations, 
and the coefficient of t' represents the sensitivity to temperature 
fluctuations. These sensitivities could in principle be calculated, 

Since all factors are either measurable quantities, physical properties, 
or known constants. In practice, however, this calculation yields values 
of the sensitivities that do not agree with those determined by actual 
calibration. This discrepancy is due to two of the supposedly known 
quantities in Equation 38, namely a and B. King's value for B, since 

it is based on his potential flow assumption, is seriously in error. The 
value of a found in tables of physical properties may also be in error 
because the value of a applicable to the wire in use depends very strongly 
On the history of that wire. In addition, contamination of the wire 


changes the value of D, the wire diameter, and the wire length & is some 





1] 


Errective length rather than the actual length. Thus it is impossible 
to accurately determine the sensitivities by direct calculation. In- 
stead the sensitivities are determined by calibration. 

Although the expressions for the sensitivities in Equation 38 do 
not provide quantitative information, they do indicate the direction of 
variation in the sensitivities as various parameters are changed. Thus, 
Beaeevelocity sensitivity should decrease with increasing velocity, but 
velocity should not affect the temperature sensitivity to first order. 
The velocity sensitivity should increase with increasing wire temperature 
(resistance). 

In order to determine the effect of parameter variation on the 
Eemperature sensitivity, further consideration of the two terms that 
comprise it is necessary. Equation 8 shows that, all other parameters 
remaining constant, an increase in qT. (i.e. a positive t') requires a 
decrease in Ee (i.e. a negative e') in order to maintain the equality. 
Thus, the overall sign of the temperature sensitivity must be negative. 


since all factors are inherently positive, this implies that 


ER a R (R -R_) nxfA 
s » W x a : (39) 
2(R -R_) 2 EaR 
w a Oo 
ER a 
Thus an increase in a causeS an increase in 
2(R_-Ro) 


mie magnitude of the temperature sensitivity, and an increase in 
R_(R _-R_)nwfa 
Ww a 


2F aR causes a decrease. 
a E 
Equation 33 shows that the quantity Rom =R ) 
wow a 


is independent of ee - T That is if the parameter q. - qT. is varied, 





sy 








a = a 
ig remains constant. Further if R. iS constant, Ee is 
- R 
RR, a) 
2 
proportional to RO Now 
E Ra R a7) “uR ee) 
O = WwW ic Om, 
Z ae R,? E R (R zie) 2 
and 
R (R -R A - - 
Ag W x, nap = FE oe a naga : 
2E aR ER R @ (41) 
O O 
° =e ~ 
aac’ - and E are independent of ae Le 
E R eek.) 


follows from Equation 40 that variation of eae has no effect on 


ER & 
O 
PR =R Co first order. Similarly it follows from 
at a? R (R_-R_) nxfA 
ww a 
feoaeton 41 that — increases linearly with 
2 EQR. 


increasing ee As previously shown, an increase in this term causes 
a decrease in temperature sensitivity. Thus an increase in eae causes 
meaecrease in temperature sensitivity. Thus it appears that velocity 
fluctuations could be measured directly by operating at very high 


temperature differences and that temperature fluctuations could be 


measured directly by operating at very low temperature differences. 





13 


If consideration is given to changing the dimensions and material 
of the wire, further effects on the sensitivities occur. Velocity 
sensitivity increases with wire length while the temperature sensitivity 
decreases. An increase in wire diameter causes an increase in velocity 
sensitivity and no change in temperature sensitivity to first order. An 
increase in temperature coefficient causes a decrease in velocity 
sensitivity and an increase in temperature sensitivity. 

Once the desired information on trends in sensitivity is recognized, 


mers convenient to resort to a more compact notation: 


le oelee ee = 1 DENG 2 
R (RU-R,)=£D Cares kk 





a » Fe 2 2 ree 
oF -& 
(42) 
and 
a ER a ; a Ky ale pan 
ti . a 
2(R -R_) 2 Ew R 
w a O 
(43) 
thus Equation 38 becomes 
J tee i S Jn 
e sou ae (44 ) 
Squaring and time averaging yields 
a > ae 
a = § oo ut tS : ao. . (45) 
u up Tae i 


The sensitivities S and Sp must be determined from calibration. 





This is discussed in the section on procedure. The quantity, e' 
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———— 


is experimentally measured. From known values of e' , Se and Sep a set 


of simultaneous equations of the form of Equation 39 can be formed. 


rr ae 





Ms 
This set can then be solved for u' , u't', and aie. This solution is 


discussed in the RESULTS section. 





lies 


APPARATUS 


F Low system 


The flow system is shown schematically in Figure 1. A Craftsman 
Model 315.16970, 3HP industrial vacuum cleaner was operated as a blower 
to supply air to the system. A three inch galvanized duct carried the 
air to the highest point of the flow system where the duct tapered to 
meet the 1.0625 inch inside diameter Type M copper tubing that comprised 
the remainder of the loop. The three-inch duct was fitted with a bleed 
line and damper. The flow rate in the loop could be varied by using the 
mamper to alter the diviSion of air between loop and bleed line. Above 
Ee Dleed line in the duct was an orifice meter for measuring air flow 
Fate in the loop. (The calibration curve for the orifice meter is 
given in Appendix B.) In the copper tube, the air passed downward 
through an 83 inch unheated velocity developing section. This was 
followed by a 58 inch heating section which was wound with nichrome 
heating wire and covered with 1.125 inches of pre-formed pipe insulation. 
Probes were inserted upward into the copper tube as shown. All measure- 
ments were made in the end plane of the heating section. 

Iron-constantan bulk thermocouples were installed at the beginning 
Or the velocity developing section and at the end of the heating section, 


ena wall thermocouples were installed at the end of the heating section. 
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Figure 1. Flow Diagram of Apparatus 
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Probes and Traversing Mechanisms 


The rapid response thermocouple used (a Heat Technology Laboratory 
TCFW 202-ChA-5) is shown in Figure 2. The response time of this thermo- 
couple has been estimated by Rodriguez-Ramirez to be 5 milli-seconds. 

The hot-wire probe used (a Thermo Systems 1210-T1.5) is shown in 
Figure 3. 

Positioning of the rapid response thermocouple was by means of the 
traversing mechanism shown in Figure 4. When in place as shown in Figure 
1 the mechanism was, of course, inverted. The ball joint permitted 
radial movement of the probe across the tube. Radial position was 
measured with the dial guage shown. 

Figure 5 shows the traversing mechanism used to position the hot- 
wire in the loop. This traversing mechanism replaced the thermocouple 
traversing mechanism at the bottom of the loop when hot-wire measure- 
ments were being made. The upper block in Figure 5 could be moved 
laterally so as to place the hot-wire probe at the desired radial position, 
while the lower plate shown in Figure 5 was bolted to the loop. A pitot 
tube and a local mean thermocouple were inserted through the auxiliary 


probe entrances. 








doe 


f= 








ey 


=? TUNGSTEN WIRE 


OOO0l5 PIAMETER) 


WIKRE SUPPORT 


Preure 3. The tlot-Wire Probe 





20 


Lial Gage 







. 
q 
sl 
le 


N 


he 






Sa Nae 
|~ " “OG Ring 
" Y Yet. Bellows 
2.563 
alr 
9 ae fy 
Y 
0625" ( 
WN 








25 —~-=: 





=o ¢ + Ses 


Fast hesponse 


thermocouple _dtog2st-— 


Figure 4. Thermocouple Traversing Mechanism 





mapure 5. 


ee 


oh@ als ihe PROBE 
AUXILIARY PROBE 


wae 





Hot-Wire Traversing Mechanism 





22 


Pustrumemtat lon 

Instrumentation for hot-wire measurements consisted of a Thermo 
Systems Type 1010 anemometer, and a Ballantine Model 320A True Root- 
Mean-Square Electronic Voltmeter. 

For rapid-response thermocouple measurements, instrumentation consist- 
ed of a Tektronix 122 Amplifier (See Appendix C for frequency response. ) 
fmcetne Ballantine Voltmeter. 

Bulk and wall thermocouple measurements were made with a Dana Model 


2850 D.C. amplifier and a Beckman Model 4OLIR digital voltmeter. 
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PROCEDURE 


Hot-wire in [sothermal Flow 
The first step was calibration. The constant C defined by Equation 
17 and appearing in Equation 18 can be determined from mean voltage and 


-2 
velocity measurements. It is seen from Equation 14 that if E 


See 

were plotted as a function of cc the result should be a straight line 
Sees slope C. 

The Calibration of the hot-wire probe was carried out with the probe 
in place in the loop and located at the tube centerline. Variation of 
u was achieved by varying the flow rate in the loop using the bleed line 
@amper. The quantity E was measured at seven or eight different flow 
rates. The value of u at the centerline corresponding to each flow rate 
was measured with a pitot tube. The data were plotted as indicated above, 
a straight line was fitted to the points, and the slope was measured. 
Calibration curves are shown in Appendix E. 

The second step was the making of fluctuation measurements. Equation 
18 gives the relationship between the rms velocity fluctuation and the 
rms voltage fluctuation. 

The flow system was operated at Reynolds Numbers of 30,700, 41,500, 
and 49,800. For each Reynolds Number, measurements were made of Vet? 
al eight points between the center and the wall of the tube corresponding 


Mme emot ©. 0.125, 0.250, 0.375, 0.500, 0.625, 0.750 and 0.875. A 
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meLot tube was used Co measure u at the same points. 
After all fluctuation measurements had been made, the calibration 
procedure was repeated to make certain that the constant C had not 


Changed during the run. 


Rapid Response Thermocouple in Non-Isothermal Flow 


The flow was adjusted so that the Reynolds Numbers in this part were 
the same as those in the previous section. The applied heating was 430 
watts in all cases. 

Ver : 

Measurement was made of e from the rapid response thermocouple 
as a function of radial position for each Reynolds Number. No calibration 
was done, Since the constant, a, in Equation 24 was assumed to be that 

ES : ; 
found in standard thermocouple tables  . This assumption is of question- 
wile validity. 

A conventional thermocouple was used to measure the centerline 
temperature, and the wall thermocouples were used to measure the wall 


temperature for each Reynolds Number. 


Hot-Wire in Non-Isothermal Flow 
The procedure of calibration and measurement outlined below repre- 
sents one of the contributions of this work. It is anticipated that 
this procedure will be used for the glycol and mercury measurements 
previously mentioned. 
As shown in the THEORETICAL BACKGROUND section, the sensitivities 


a and Sp depend on two parameters: u and yO 3 [ivess anv calibration 


* 


it 1S necessary to determine the functional dependence on each parameter. 


Calibration for temperature sensitivity was also carried out with 
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mieeprobe in place and at the tube centerline. Velocity variation was 
again achieved by use of the bleeder line damper, and u was measured with 
a pitot tube. In this way five different values of u were used, and 

for each value of u thirteen values of ete were used. 


These data were plotted as E versus u for constant values of z =n 


and as E versus z as for constant values of u. These curves are shown 
Maeeppendix G. The sensitivities S and Sp are the derivatives of these 
curves. These derivatives were obtained from analytical differentiation 
of parabolas fitted to the data by the method of least squares. The 
selection of parabolas is arbitrary, but the relatively slow and smooth 
variation of the data makes a parabola a satisfactory approximation. 

Note that ue must be computed from measurements made of R-Ro: 
The resistance-temperature characteristic of the hot-wire probe was 
experimentally determined to be as shown in Appendix D. 

Once S and Sp were known as functions of u and i To fluctuation 
measurements were made. The rms voltage fluctuation was measured at 
r/R = O and r/R = .5 for each flow rate and temperature difference con- 
Sidered in calibration. 


When all fluctuation measurements had been completed, the calibration 


was repeated to determine the drift in oF and Sa 
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RESULTS 


Hot-Wire in Lsothermal Flow 

Preare 6 Shows the axial turbulent intensity at the centerline of 
the pipe compared with the work of Sandborn. 

Figure 7 shows axial turbulent intensity as a function of radial 
position for three different Reynolds Numbers compared with the work of 
Reteer 

The agreement of present results with those of SECC and 
fre” is rather good. The data used by Sandborn to obtain his 
correlation show a scatter of + 10 percent. The results of the present 
work are well within this range. 

Rapid Response Thermocouple in 
Non-Isothermal Flow 

Figure 8 shows the temperature fluctuation intensity as a function 
Semradial position. eee = hamirena. and also Tanimoto and Hanratty 
have published results of this nature. The present results are for a 
range of Reynolds Numbers higher than the range considered by Rodriguez- 
Ramirez; however, a preliminary run was made in which the flow rate and 
heat input were the same as those in one of Rodriguez-Ramirez' 
runs. The results of this run deviated from those of Rodriguez-Ramirez 
by less than 5 percent. The temperature gradients used by Tanimoto and 
Hanratty were so much larger than those used in the present work that 


comparison is not possible. 





Figure 6. 
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Hot-Wire in Non-Isothermal Flow 

Figures 9 and 10 show the dependence of the sensitivities S and Sp 
Saeene parameters u and eT The velocity sensitivity, Se increases 
with increasing temperature difference and decreases with increasing 
velocity as predicted from Equation 38. The temperature sensitivity 
decreases with increasing temperature difference as predicted from Equation 
38. Note that there is a secondary dependence on velocity not predicted 
by Equation 38. This dependence is due to the velocity dependence of the 
Serective wire length. 

The investigation of various calculating schemes including the 
selection of the best one for analyzing the non-isothermal data was an 


mmportant part of the present work. As pointed out in the section on 


theoretical development, Equation 45 gives the relationship between the 


—_ —S ee 





2 
unknown quantities u' , u't' and t' and the known quantities Se Sip 


and Eo . Since there are three unknowns, at least three equations of the 
form of Equation 39 are needed. These can be obtained by operating the 
hot-wire at least three values of ie (and hence of Ie For each 


value of 1 elie the values of Ss. and Se appropriate to the mean velocity, 


u, of the flow are read from Figures 9 and 10, and pg is measured. Thus, 
for each value of ee an equation is formed. 

Once the set of equations is obtained, a number of methods of solution 
are available. The most obvious one is the simultaneous solution of 
three such equations by a method such as Gaussian neaue lente However, 
the numerical calculations involve the subtraction of large numbers whose 
difference is small, which has the effect of greatly magnifying the errors 
in the data. Since the errors in Sy and So were already large, the result 


Sm ee ee 


em 2 
mimeorrors in the calculated values of u' , u‘'t', and t' of the 
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order of 500 percent. 
; : LES 
The Kovasznay Fluctuation Diagram is a graphical method of solving 


Z 
Bre set of equations in which equation 39 is divided by D1 to give: 





72 5 rns S ——5 
S_—=(—) ao) ee (46) 
S ; am oT 
i 
et? S. 
Maen the quantity go is plotted versus a Note that analytically 
S T 


iT 
Mae relationship between the two is described by a parabola. By the 
method of least squares a parabola is fitted to the data. The computer 


program used to do this fitting is listed in Appendix F. The intercept 








12 S 
Bietnis parabola with the a axis is rte eee Ss hOpemas = =0is 
Se2u t' , and the curvature of the parabola is u'  . This method was 
S 

Msed on the present data, but the need to extrapolate the data from oe of 

5 5 T 
Some non-zero value to  G = QO and then determine the slope at — = 0 

T T 
resulted in a larce uncertainty in u't’ . This same problem was encountered 


2 ’ 
by Arya and Plate , who suggested that if an independent measurement of 


2 
' could be obtained to locate the intercept of the curve accurately, 








ae 2 
then reliable values of the remaining two unknowns u't!' and u'” could 


be determined by the Kovaszgay method. This procedure was followed with 





Phe present data. The value of ue used was the value found using the 
mma response thermocouple. The results were still rather uncertain, 
but the uncertainty was less than for any of the other methods. (The 
Kovasznay Fluctuation Diagrams for representative cases are shown in 


Appendix H.) 


I 


NO 


The numerical values of ¥ u' Ju ¢. presented in Figures 11 and 12 


deviate seriously from the isothermal resulls wich which it is thought 


es 
they should agree. (See, for instance Corrsin and Uberoi ') The trends 
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of the results with Reynolds Number and radial position are in accordance 
with expectations. 

The calculated values of lege el ane presented im Figure 132 
Turbulence at the centerline in pipe flow for the Reynolds numbers used 
in the present work should be isotropic. Thus u't' and aay ie (T-T) 
should both be zero at the centerline. The centerline results in 
Figure 13 do not exhibit this behavior. However the error in these re- 
sults is roughly estimated in Appendix I as 100%. 

Away from the centerline the turbulence should no longer be iso- 
PEOpic , and u't! should no longer be zero. With heating from the wall 
as in the present case, u't', and hence u't'/u, (T -T.) should be negative. 
This trend is followed by the present results in that u't'/u (tT -T,) at 
r/R = .5 has a larger magnitude than at r/R = O and in that it has a 
Mesative sign. The quantity tiene a cannot continue to increase 
indefinitely as Reynolds Number increases, for this would mean no dis- 
Sipation was taking place. The increasing trend in the results shown in 
Breure 13 for r/R = .5 is therefore not anticipated, but it is possible 
So long as the trend reverses itself at higher Reynolds Numbers. Arya 
and iste” found that in the boundary layer growing on the wall of a wind 
tunnel u't'/u, (T.-T, )remained nearly constant with increasing Reynolds 
Number. Again with an estimated error of 100 percent little significance 


can be attached to present results. 
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DISCUSSION 





The numerical results for V ut /a, and cena 
are in general, too large. In view of the way in which these intensities 
are obtained from the Kovasznay Fluctuation Diagram, this suggests that 
Sy iemtoo Small. The temperature sensitivity Sap affects both coordinates 
in the same way, and therefore, errors in Sp would be self-compensating. 


S 


: u ; 
mowever, errors in 3 affect only the = coordinate. Too small a value 
T 


of aa would cause too large a slope and curvature in the fitted parabola 
t 2 1 t 

and too large a value of u shoal El [tse 

Because of the long time required for the system to come to thermal 
equilibrium after each change in flow rate, the data used to compute 5$ 

u 
mocea Siven value of Teel. were taken over a period of approximately 
eight hours. It is reasonable to expect an increase in the thickness of 
the layer of contamination on the wire during this period since the wire 
meas Continuously exposed to an inadequately filtered flow. 
3 . ; 

Kronauer has shown that dirt accumulations on fine wires at low 
flow rates increase the heat transfer rate because the heat transfer area 
[eereases more rapidly than the heat transfer resistance as dirt 


accumulates. Thus, a dirty wire would require a higher value of E than 


eeclean one to maintain it at a given temperature. 





ae, 


The wire used for the present measurements was examined under a 
microscope, and a significant amount of dirt was observed on it. 

manene present calibration, E was measured at the highest calibration 
velocity first, and the velocity was decreased step by step until the 
lowest value was reached. As anticipated, E decreased with decreasing 
velocity, but the increasing layer of contamination apparently worked 
counter to this trend. Thus, in the calibration procedure used the rate 
of decrease of E was slower than if contamination had not been a prob lems 


Maes result was a value of 3 ENaAte was too cma ll. 





= 
Larger values of ee would result in larger values of = at the same 
2 ve 
value of e et on the Kovasznay Fluctuation Diagram and thus, a 
2 
> 
T 


smaller slope and curvature of the fitted parabola. The calculated values 
V 2 Pea cat . ; 
of u and u't' would then be smaller and in better agreement with the 
meochermal results. 
. - 12 sage 
HE was not possible to determine the errors in u and ut 
Caused by the drift in E noted above because it was not possible to 


measure what E would have been without contamination. Thus it was 


mecessary to use a different method to estimate these errors. The method 





Meed is outlined in Appendix I, and the error estimates are 36 percent 
in Vue and 100 percent in u't'. As noted in Appendix I, these figures 
Mommot represent true confidence limits on the results. They are only 
rough estimates of the errors present. 

The calibration curves for before and after as shown in Figures 24 
Phrough 26 show such a large drift that all data taken in between should 


be discarded. This was not done, since a number of attempts to obtain 


 t=* — 4 
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better results failed. Note then, that all fluctuation intensities ob- 
tained from the present non-isothermal hot-wire measurements are of very 
doubtful reliability. 

Due to inexperience in this area, a number of deviations from good 
Seperimental practice were present. First, although the wire was 
examined under a microscope after all measurements had been made, no 
examination was made before or during experimental runs. Thus the 
hypothesis of an increasing dirt layer cannot be verified. Second, no 
attempt was made to clean the wire during the course of experimental 
runs. Because of these errors it is not possible to analyze in detail 
Memeo correct for the drift that occurred. 

In future work, provision should be made for careful filtration of 
the fluid and a means of rapid calibration to minimize drift. Two cali- 
brations, one with increasing velocity and one with decreasing velocity, 
should be made to check for contamination. If it is still not possible 
Memreamee drift to an acceptable level, provision should be made for the 
complete withdrawal of the probe from the flow. 

A problem encountered with the present apparatus was difficulty in 
maintaining flow rate and heat input at desired levels. Long term 
memperature variations of two or three degrees persisted even after the 
system had come to "equilibrium". This was apparently caused by cycling 
of the room temperature. 

The nature of the problems encountered in the present work serve as 


feeuide in planning additional steps in the overall research program. 
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SUMMARY AND CONCLUSIONS 


This work has developed procedures for the calibration of a hot- 
wire in non-isothermal flow. Further, the various calculation schemes 
available for analyzing non-isothermal flow fluctuation data have been 
evaluated and the Arya and plea” modification of the Kovasznay Fluctuation 
eran | method has been found to be the most satisfactory. 

Apparatus limitations and deviations from good experimental practice 
meeciuded the production of reliable numerical results. 

The findings of this investigation are an important preliminary step 
in a comprehensive research program of non-isothermal flow measurements 


with hot-wire and hot-film anemometers. 
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APPENDIX A 


TABLES OF )XPERDMENTAL DATA 


Table 1] 
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Table 3 


Calibration for Isothermal Re=49,800 Run 
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Table 6 


Data for Isothermal Re=30,/00 Runs 


Radial Guo hPa i Pot“ tube - Fa 

Position Crm so PG ./797 © paid) E (Wolts) v e'° (Millivolts) 
0.0000 1.480 Zouk hee 2 
On0625 Ic460 2eol V2.9 
e250 1.400 2.80 eal 
O-18/5 eee 6, 2.79 oe 2 
2 500 1.240 Za Zao 
Oe 3.125 e130 aa is) 24.5 
3 / 5.0 PAOZ0 2274 28.6 
024375 0.850 2.71 B25 

Table / 
Data for Isothermal Re=41,500 Run 

Radial BE Pittrot. Tube = 

Position Ce (ties F679 7a iid) E (Volts) v e' (Millivolts) 
O20000 L130 295 254 
Oe0625 famed Oe. 2.95 Pock 
0.1250 22650 ZOOS ego 
Ole 75 2.500 2.94 19.0 
O25 00 ZS 30 pi is 276 
Oes25 Vales, (0. 2790 On 2 
O37 50 1.960 Vie eh) 29.4 
074375 le Ga0 22Go) bee 2 
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Tables. 
Thermocouple Data for Re=30,/700 Run 


Radial 
Position 
fin ) 


Pea! 2 
(Milli- 
Meme 20.0 21.0 23.5 28.5 30.5 34.5 38.5 42.5 45.0 47.5 8.0 


OOO Or OommO Om 0, 15) OF 20e80725" 0230" 0.35 0.40. 0.45 NOISE 


Table 9 
Thermocouple Data for Re=41,500 Run 


Radial 
Position 
cin ) (ae OeOsmr Oo TO Oe la 0.20 0.25 0230 0:35 -0.40' 0.45 NOISE 


7 at? 
(Milli 
meee.) 15.7 18.3 21.5 25.0 28.5 32.0 35.5 38.0 39.5 6.0 





Table 10 
Thermocouple Data for Re=49,800 Run 
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y eZ 
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Non-Isothermal Hot-Wire’ Before Calibration Data 
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Table Ve 


Cold Resistances in'8efore'Calibration 


P(ina.) = R eo ig (OHM 
£220 7:42 
io 20 7.46 
io30 7.49 
05932 rey: 


0.480 re 18. 
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Table 13 


Temperature Difference Between Wall and Centerline 


_ ain) © | [R=] ; T r/R=0 a 
22420 46 
1.620 Dt 
1.330 54 
ee 61 
.480 74 


*Orifice Meter Manometer Fluid SPG 2.95 
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Table 14 


Non-Isothermal Hot-Wire RMS ier oe 


Fluctuation Data = 


7 AP (in)*-+ 
“FE (Ohm) ese i520 ies sOF 0.932- 0.480 
{ 






OO Mil lae iad 





volts 
ie 2 0 24.8 EC 2324 SG, OTe 
1.40 Zak] iLowee ie ies SOS: BO 2 
PO O 24.5 201. 0 299 ooo 33:53 
1.80 HEISE ee) 2 font ee sys 
2.00 Zo eS ee. 26729 28.9 33.0 
20 Cane CAS a? 2Or? 20 <5 Bed 
2.40 Zee il Zo.) 2604 ate 32.0 
2.00 24.1 29 2642 tae eee aia: 
230 LESS Coe] 25.9 ve GE 3150 
Bb 00 24.0 Cas oo wets tad OP 5056 
B20 se Dees, 24.5 257.5 ara’) 30-0 
Be +0 259 eo 2D PAY Zee 
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Table 15 


Non-Isothermal Hot-Wire RMS Voltage 
Pilvetwation Data tor r7kR=.5 


AP(in)* >» 
R -R (Ohm) 
oa 


t 


egy ee N78, ieee 0.932. “09480 





Milli- 


| Oe) G27 2s Volts 45.20 46.8 a0c Disa 
ie 20 aaa 44.8 46.8 20.0 35.3 
re 0 42.6 44.6 46.2 49.2 see 
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Lite) Pa 44.0 45.6 48.0 Do0 
2.00 “2.7 43.8 ae re | 47.5 54.2 
ey, 20 253 43.5 44.9 47.0 pares. 
2.40 423 43.3 44.8 46.2 Dee 
= 00 ees aa 44.5 46.0 Blew 
OO 42.8 Le 44.5 Ls 3 dO ete 
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B20 42.6 43.0 44.3 oral 49.5 
3.40 e216 43.0 44.0 44.8 49.2 


*™ Orifice Meter Manometer Fluid SPG 2.95 
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Table 16 


Non-Lsothermal Hot-Wire'After"Calibration Data 





AP Ca) 
oe, (Ohm) a a Les 20 be330 Oey 0.480 
1.00 goat, (PSE 1,945 VAG. 1560 Le Be 
1.20 2.140 BOO. <2:.07.0 see hes 12670 
PO 22235 ag he Ze lo0 JNA G, 2.000 
i 00 22305 PBA, 2285 22220) eles 
ie 3 0) 2.470 2.020) 2235) Ze) eS 
2.00 22909 ae a 2.480 2 aS 7 ger ae Ns) 
20 2.695 22000) 22565 2.495 ares Wee. 
24.0 2.740 25005 Ze Jamal Die) 
2,00 27620 2 Oo Pe 20S 2645 Zoe 
2 OU 2.900 2.840 2.795 aan Boe W Lay Fs 
3.00 SSE, 7p 2) bs) Peco)? 2755 ya ops 
20 3070 2.980 22930 8 €2..045 Cals 
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Meorifice Meter Manometer Fluid SPG 2.95 





Table 1/7 





Cold Resistances in'After''Calibration 


Pein.) 7 R old (OHM) 
2b 20 ra! 
1.620 7.24 
1330 P26 
Or 732 Tiere 


0.480 ae) 








AP PEND TA 








a4 


APPENDIX B 


ORTPICE METERS CALIBRAT LON 


The orifice meter used was calibrated in place against a pitot tube 
located at the centerline of the test section. von Karman's universal 

: fees! | 
Meeeecity distribution was used to calculate flow rates from the pitot 
tube velocity measurements. The calibration curve is shown in Figure 14. 
Meeis tound as anticipated that a plot of flow rate versus the square 


mor of the pressure drop is a straight line. 
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APPENDIX C 





APPENDIX C 





AMPEIP IER FREQUENCY RESPONSE 
The frequency response of the Tektronix amplifier is shown on 


the followin & Peet. 


DY 
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APPENDIX D 


RES tot ANCE -TEMPERATURE CURVE FOR THE HOT-WLRE 


Figure 16 shows the relationship between resistance and temperature 


Mer tie hot wire used in the non-isothermal flow measurements. The 


@ashned curve is the approximate relationship suggested by Kronauer. 
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APPEND EXE 
CALIBRATION OF HOT-WIRE IN [SOTHERMAL FLOW 


The isothermal flow data were taken on two diftterent days. The set 
of data for a Reynolds Number of 49,800 was taken on the first day. The 
@eiypration curve for these data is shown in Figure 1/7. The constant C 
mamecmation 18, which is the slope of the line in Figure 17 is ier 10 

The data at Reynolds Numbers of 30,700 and 41,500 were taken on the 
meeond day. The calibration curve for these data is shown in Figure 18. 
mie constant C tor these data is eeerion 

immemyalue of C calculated from Equation 17 is 0.84 x aii. The 
discrepancy between this value and that obtained from calibration is due 


Morwire history effects in ® and R_ and the error inherent in King's 
. O 


potential flow value for B. 
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APPEND i 


LISTING OF CURVE FITTING PROGRAM 
The computer subroutine used to perform least square curve fitting 


calculations is listed on the following pages. 
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APPENDIX G 


CALIBRATION OF HOT-WLRE IN NON-ISOTHERMAL FLOW 


Figures 19 through 23 represent the calibration data used to compute 
the results presented for non-isothermal flow. The curves drawn through 
the data represent the least squares fit of parabolas to the data. These 
parabolas were differentiated analytically to obtain the values of and 
Sr presented in Figures 9 and 10. 

Figures 24 through 26 show calibration curves for representative 
values of u and ee - i before and after fluctuation measurements were 
made. Table 18 gives a complete tabulation of S and 5, as functions of 


T 


u and sm = v before and after the fluctuation measurements. 
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Figure 19. Non-Isothermal Calibration of Hot-Wire I 
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Figure 20. Non-Isothermal Calibration of Hot-Wire II 
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KOVASZNAY FLUCTUATION DIAGRAMS 


Figures 27 through 29 show the Kovasznay Fluctuation Diagrams for 


B@eeesrepresentative cases. 


ime parabolas that fit these data are: 








Preure Parabola 
ae ae 5 Coe 2 
27 ( ~t'°) = 2.64 +16.9(—) + 19.1(—) 
2 S S 
S T db 
dh 
ae ae Sy Sh Z 
28 ( =~ - t' ) = 0.43 +18.9(—) + 16.2(—) 
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5) T i) 
T 
et? 9 S S 2 
29 (== - t'*) = 0.21 + 42.6(—) + 47.3(5 
2 S 5 
5 T il 
ih 
SF 
a 1S then Lound by dividing the coefficient of ae by <2. 
a T 7 2 
and  ¥ u' toweme scdudre sw Oolwol ihercocertieient of  <( a Pati 
Tt 


results from the parabolas given above are listed in Table 19. 
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Figure 27. Kovasznay Fluctuation Diagram I 
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Figure 28. Kovasznay Fluctuation Diagram II 
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Figure 29, Kovasznay Fluctuation Diagram III 
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APPENDIX I 
ERROR ANALYSIS 


The results of the isothermal flow measurements with the hot-wire 
deviate from the results of See ae and waeeee by less than 10 
Penreent. 

The results of a preliminary non-isothermal run duplicating the 
conditions of Rodriguez-Ramirez and using the rapid response thermocouple 
agreed with his results within 5 percent. 

Since the errors in the results of the non-isothermal flow measure- 
ments with the hot-wire are quite large, a detailed estimate of their 
Magnitude is needed. 


t 
u 


Do 2s ——__ _ 
It is assumed that the errors in V /u. and ult'/uy (T .-T.) 


& 
are primarily the result of errors in oy and Sep In order to obtain an 
estimate of the errors in 5 and So an analysis was made of the dis- 


crepancies between successive calibrations as tabulated in Table 18. 
Note that the actual errors in Sh and Sp come from a continuous 
deitt in E Ctiainesasciven Calibralion. as outlined in the DISCUSSION 
Section. Since it iS not possible to meaSure this drift, recourse was 
made to analyzing by this method. Thus the error estimates obtained 
BeloweGo Nek represent.contidence limits on the results but only rough 


SypGoximationc ste the contidence limits. In Table 18 , the largest 
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discrepancies occur for u = 137 ft/sec. The percentage error, based on 
the before value, between the before and after values, was determined 
for each entry in the 137 ft/sec column. All the S.. errors were summed 
and divided by the number of entries and likewise for We Thus, average 
errors were determined to be 19 percent in S and 14 percent in Sip 

The propagation of these errors into the coordinates of the Kovasznay 


Fluctuation Diagram is as follows: 








2 
6 Ae = ae 
> z as) 
iT. east = 28 percent. 
= 9 
eat 2 : 2 li 
Z 
S 
18 
5 
ACU) 
a 
oT - a + aA = 33 percent. 
2 i “T 
@) 
ah 


These errors were assumed to apply to all values of u and r/R. 
5 
The value of a for each point of the KovasZnay Fluctuation 
T 
Diagram for u = 66 ft/sec r/R = 0 was increased by 33 percent, and the 


2 
value of ie =a was decreased by 28 percent. A parabola was fitted 
Z 
oT 


Y. ere 
to the resulting points and u' and u't! determined. The new values 
differed from the original values by 72 percent and 100 percent respec- 


pivelv. Since 


a 
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hole 
Cc 
5 
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2 
the ericor WAN is 36 percent. This 36 percent error was used as 





Z 
Ene limit of error for all values of NITE pod , and 100 percent was used 


ee ; 
foemeall values of u't /u, 7, qT) 
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